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Abstract 



Bottom quark production in pp collisions at s/s = 1.8 TeV is studied with 5 
pb _1 of data collected in 1995 by the D0 detector at the Fermilab Tevatron 
Collider. The differential production cross section for b jets in the central 
rapidity region (\y b \ < 1) as a function of jet transverse energy is extracted 
from a muon-tagged jet sample. Within experimental and theoretical un- 
certainties, D0 results are found to be higher than, but compatible with, 
next-to-leading-order QCD predictions. 
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Measurements of the bottom-quark production cross section at pp colliders provide an 
important quantitative test of quantum chromodynamics (QCD). The mass of the b quark is 
considered large enough {m b ^> Aq CD ) to justify perturbative expansions in the strong cou- 
pling constant a s . Consequently, data on b quark production are expected to be adequately 
described by calculations at next-to-leading order (NLO) in a s [0 0]. 

Past measurements of inclusive b quark production in the central rapidity region, at 
center-of-mass energies yfs = 0.63 TeV |||| and ^Js = 1.8 TeV indicate a general 

agreement in shape with the calculated transverse momentum (p?) spectrum, but are sys- 
tematically higher than the NLO QCD predictions [TJ-[J by roughly a factor of 2.5. More 
recently, a measurement using muons at high rapidity |[0) (2.4 < \y^\ < 3.2) to tag b 
quarks, indicates an even larger excess of observed cross section over theory. The measured 



differential production cross section for B mesons [|TTJ is similarly higher than the NLO 
QCD prediction. Calculations of higher order corrections [|l|,|l3|] have shown that additional 
enhancements to the cross section beyond NLO are likely. However, the expected enhance- 
ments fall short of accounting for the observed discrepancy between theory and data. This 
longstanding mismatch has motivated continuous theoretical and experimental effort dedi- 
cated to reduce uncertainties in general, and to broaden the scope of observable quantities 
Ml. 



Previous studies of b quark production by the D0 collaboration have exploited the kine- 
matic relationship between b quarks and daughter (semileptonic decay) muons to extract 
integrated b quark production rates as a function of p h T threshold J7j, and have examined 
azimuthal correlations between the b and b in pair production || . 

The present study is a complementary measurement of b production, based primarily on 
calorimetry, with the main focus being on b jets rather than b quarks, b jets are defined 
as hadronic jets carrying b flavor. As opposed to quarks, jets are directly observable and 
therefore reduce model dependence when comparing experimental data with theory. This 
measurement is in direct correspondence with a NLO QCD calculation [I5| that highlights 
the advantages of considering b jets rather than open b quarks. For instance, large logarithms 
that appear at all orders in the open quark calculation (due to hard collinear gluons) are 
avoided when all fragmentation modes are integrated. 

The differential production cross section of b jets as a function of the jet transverse 
energy (E T ) has been measured using data collected during 1995 with the DO detector 
at the Fermilab Tevatron Collider. The data correspond to an integrated luminosity of 
(5.2 ± 0.3) pb _1 . The analysis exploits the semileptonic decays of b hadrons which result 
in a muon associated with a jet. An inclusive sample of muon-tagged jets is selected from 
a trigger requiring a jet and a muon, and the b jet component of the inclusive sample is 
extracted, based on the properties expected of the associated muon and jet. 

The DO detector and trigger system are described elsewhere Jet detection utilizes 
primarily the uranium-liquid argon calorimeters, which have full coverage for pseudorapidity 
|j~7| \r]\ < 4, and are segmented into towers of At] x A<ft = 0.1 x 0.1, where <ft is the azimuthal 



angle. The relative energy resolution for jets is approximately 80% E(GeV). The central 
muon system consists of three layers of proportional drift tubes and a magnetized iron toroid 
located between the first two layers. The muon detectors provide a measurement of the muon 
momentum with a resolution parameterized by S(l/p) = 0.18(p — 2)/p 2 © 0.003, with p in 
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GeV/c. The total thickness of the calorimeter and toroid in the central region varies from 
13 to 15 interaction lengths, which reduces the hadronic punchthrough for the muon system 
to less than 0.5% of muons from all sources. 

Initial event selection used a trigger requiring (i) Ex > 10 GeV in at least one calorimeter 
trigger tile of Arj x A0 = 0.8 x 1.6, (ii) at least one muon candidate with px > 3 GeV/c, 
and (iii) a single pp interaction per beam crossing as signaled by the trigger scintillation 
hodoscopes. 

Jet candidates were reconstructed offline with an iterative fixed-cone algorithm (cone 
radius of 0.7 in rj-cf) space) and then subjected to quality selection criteria to eliminate 
background from isolated noisy calorimeter cells and accelerator beam losses which may 
mimic jets. Accepted jets were required to have Ex > 25 GeV, \r)\ < 0.6, and an associated 
muon within the reconstruction cone. 

Offline muon selection required the triggered track to originate from the reconstructed 
event vertex, with p T > 6 GeV/c and [rj^l < 0.8. Muon candidates in the azimuth region 
80° < ^ < 110° were excluded due to poor chamber efficiencies near where the Main Ring 
beam pipe passes through D0. 

Background events from cosmic ray muons were eliminated by requiring a 10 ns wide 
time coincidence of muon tracks with beam crossings. About 30,000 events were selected, 
of which less than 0.5% have either a second muon-tagged jet or a double-tagged one. This 
is in agreement with expectations from muon cuts and acceptance, branching fractions, and 
flavor content of the sample. 

The differential b jet cross section as a function of jet Ex is extracted from the inclusive 
tagged jet sample as 

da U(E T )F b (E T ) AN» J 



dE b T ' et 2BC int e flJ (px, Ex)A(Ex) AE T 

where AN^ -7 /AEx represents the inclusive tagged jet counts per transverse energy 
bin-width, Fb(Ex) is the fraction of tagged jets containing muons from b quark decays, 
e Mj (py, Ex) is the event detection efficiency of a /i+jet pair, U(Ex) corrects for spectrum 
smearing from jet energy resolution, A(E T ) corrects for the muon tagging acceptance, £ int 
is the integrated luminosity of the data sample, and B = 0.108 ± 0.005 is the branching 



fraction for inclusive decays b — > /i + X [|T^]. Tagging muons of both charges are counted, 
and the number of events is divided by two. The cross section is determined for the (b + b)/2 
combination. In Eq. |0.1| and throughout this paper, Ex represents the calorimeter-only com- 
ponent of the jet transverse energy: it excludes the tagging muon and associated neutrino. 
Once the correction for the lepton energy is included (see below), then E^ measures the 
complete b jet. 

The trigger and offline reconstruction efficiencies were determined from Monte Carlo 
events, and crosschecked with appropriate control samples of collider data. Monte Carlo 



events were generated with ISAJET |TS[, followed by a GEANT |20j simulation of the D0 
detector response, trigger simulation, and reconstruction. The overall muon geometric ac- 
ceptance multiplied by efficiency increases from 35% at p T = 6 GeV/c to a plateau of 45% 
above 10 GeV/c. The overall jet efficiency increases from 85% at Ex = 25 GeV to a plateau 
of 97% above 45 GeV. 
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FIG. 1. Fraction of muon-tagged jets which come from b decays. Et is the calorimeter-only 
component (see text) of the b jet transverse energy. 



The transverse energy of each jet was corrected for the underlying event, additional 
interactions, noise from uranium decay, the fraction of particle energy showered outside the 
jet cone, detector uniformity, and detector hadronic response [^TJ. The steeply falling Et 
spectrum is distorted by jet energy resolution, and was corrected as described in . 

In addition to b production, muon-tagged jets can also arise from semileptonic decays 
of c quarks, and in-flight decays of it or K mesons. Muons from Drell-Yan or on-shell 
W/Z boson production are not expected to have associated jets, and Monte Carlo estimates 
confirm a negligible contribution. The background from light flavors is suppressed using the 
transverse momentum of the muon relative to the associated jet axis (p^ 1 ) as discriminator. 
The higher mass of the decaying b quark generates a rather different p™ 1 distribution than 
obtained from quark decays of lighter flavor. 

Individual p™ 1 distributions for b, c, and tt/K decays were modeled using the ISAJET (19| 
Monte Carlo, with individual samples of each flavor processed through complete detector, 
trigger, and offline simulation. Due to the similarity of p^ 1 for c and tt/K sources, within 
resolution these distributions were indistinguishable. The b jet signal was extracted on a 
statistical basis, through maximum likelihood fits to the observed p^ 1 distribution. The 
normalizations of Monte Carlo templates representing light (c and tt/K) and heavy (b) 
quark decay patterns were floated such as to fit the observed p™ 1 spectra in four ranges of 
transverse energy. The extracted fraction of b jets in each Et range is shown in Fig. [I]. The 
systematic uncertainties were estimated by varying within errors the input distributions to 
the maximum likelihood fits. 

The parameterization of the overall fraction of b jets in the inclusive sample as a function 
of jet Et was obtained from a two-parameter fit (% 2 = 0.96) to the four bins in Et, and 
is shown in Fig. [I], together with the uncertainty band (6% relative uncertainty at 30 GeV, 
39% at 100 GeV). The form (a + (3/E T ) was chosen after Monte Carlo fit trials. There are 
partial bin-to-bin correlations in the uncertainty. The extracted b fraction is the dominant 
source of uncertainty in the cross section. The upper reach in b jet Et is limited by the 
deteriorating discrimination power of p^ 1 with increasing Et- 

Corrections for acceptance loss stem from (i) the tagging threshold for muon pt and 
range of pseudorapidity, (ii) the muon-jet association (tag) criterion, (iii) the undetected 
lepton energy of the jet, which is added (statistically) to the hadronic E T registered in the 
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FIG. 2. Differential cross section for b jet production. 



calorimeter, and (iv) the restricted pseudorapidity range for jets and muons (theory fl"5|l 
considers b jets with < 1). 

The acceptance correction is extracted using the ISAJET |19[ Monte Carlo simulation, 
and is defined as the ratio of b jets satisfying the analysis conditions to b jets satisfying the 
cross section definition. The transition from Et to E^ et is thus included in the correction. 
The simulated Et and rj distributions are in agreement with the data. The overall correction 
increases from 6% at E T = 30 GeV to 10% at 100 GeV. It is independent of assumed event 
production rates, and is affected primarily by models for fragmentation and decay. Decays 
are based on the QQ decay package [^] from the CLEO experiment. ISAJET uses the 
Peterson fragmentation model [p^j , with a parameter e = 0.006 that is varied by ±50% to 



estimate uncertainty. It is observed that the impact of fragmentation on the acceptance 
is mainly due to the migration of tagging muons across the minimum-p^ boundary. This 
uncertainty propagates to the cross section as a 9% effect independent of jet Et- 

The resulting b jet cross section for |?7| < 1 is shown in Fig. 0, together with the NLO 
QCD prediction |jl5| . Inputs to this calculation are the renormalization scale (chosen equal 
to the factorization scale), the b quark pole mass (4.75 GeV/c 2 ), the parton distribution 
function (MRSA' p5|), and the parton clustering algorithm. The uncertainty in theory is 
dominated by the QCD renormalization scale dependence, with the central value in Fig. |2| 
chosen as /i = [E\ + M 2 ) 1 / 2 , and varied between 2fi (lower curve) and /io/2. The fact 
that there is such a strong dependence on scale suggests that important higher-order terms 
in the calculation are still missing. 

The overall systematic uncertainty in the cross section has contributions from the in- 
tegrated luminosity (5%), trigger (3%) and offline selection (4%) efficiencies, jet Et scale 
(15%-8%) and resolution effects (7%-3%), tagging acceptance (9%-12%), and the extracted 
b fraction of jets (6%-39%). Bin-to-bin errors are fully correlated for all sources of uncer- 
tainty, except for partial correlations in the b fraction of jets. Cross section values and overall 
uncertainties, as plotted in Fig. ^ are listed in Table |I[ 

Figure |2| displays the same general pattern of past b production measurements PHIIl, 
with data lying above the central values of the prediction, but comparatively less so in the 
present case, where general agreement between measurement and the upper band of the 
theoretical uncertainty is observed. 

To connect the present measurement with previous findings, and in particular the appar- 
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TABLE I. Differential cross section for b jet production. 



rpbjct 

Cjrp 


j_ / J rpbjet 

an J arjj. 




U ncert aint ies ( % ) 


{(jev ) 


^no/ijeV J 


stat. 


oyst. 


27.3 


14.4 


1.2 


20.6 


32.3 


6.57 


1.6 


20.1 


37.4 


3.22 


2.1 


21.2 


42.4 


1.65 


2.7 


23.2 


47.4 


0.932 


3.4 


25.6 


52.4 


0.530 


4.2 


28.0 


57.4 


0.292 


5.4 


30.2 


64.7 


0.174 


4.6 


33.2 


74.7 


0.0678 


6.8 


36.6 


84.7 


0.0364 


8.7 


39.3 


94.6 


0.0176 


11.7 


41.5 



ent difference in normalization with respect to theory, the present data sample (muon-tagged 
jets) can be re-analyzed from a different perspective. Instead of focusing on the differential 
b jets cross section, the same data can be used to reproduce the integrated b quark cross 
section as a function of minimum p\. Similar measurements are described in Refs. [| [7J 
and ||. The kinematic relationship between daughter muons (pj. spectrum) and parent b 
quarks (p b T spectrum) is used to extract the b quark production cross section, integrated 
from p™ in to infinity, and over the rapidity range \y b \ < 1. Here p™ m is defined as the 
allowed p\ threshold for a given pj> distribution. Since the b quark signal fraction increases 
with muon p T (decreases with jet E T ), b tagging based on p™ 1 as a function of p^ is more 
precise than as a function of jet E?- The model dependence, introduced in the conversion 
of the muon p T spectrum to the integrated b quark spectrum, dominates the uncertainty in 
this measurement. 

The results of this re-analysis are shown in Fig. |3] (triangles) along with the theoretical 
expectation [1]] and the previous D0 |J results. Consistency among all measurements is 
evident. The present analysis extends significantly the p\ reach of the previous measure- 
ments. This is due to the requirement of jets in the trigger and offline selection, as opposed 
to requiring just muon triggers. The jet requirement naturally forces the sample into a 
higher region ofp^. Fi gure |3] shows that the agreement between data and theory improves 
somewhat at higher p\. 

In conclusion, two distinct measurements of b quark production in a previously unex- 
plored region of larger transverse momenta have been presented, and are found to be com- 
patible with the NLO QCD calculation for heavy-flavor production. Agreement between 
data and theory, unsatisfactory for previous measurements, is now observed to improve at 
higher transverse momenta. 
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this work, and acknowledge support from the Department of Energy and National Sci- 
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NLO QCD, MRSR2 

Theoretical Uncertainty 
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